Optimal parametrizations of adiabatic paths 



J. E. Avron/ M. Fraas/ G.M. Graf,^ and P. Grcch^ 

^Department of Physics, Technion, 32000 Haifa, Israel 
^ Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 
(nDated: May 19, 2010) 

We show that, given a path of Hamiltonians and a fixed time to complete it, there are many 
parametrizations, i.e. time-tables for running along the path, such that the ground state of the initial 
Hamiltonian is mapped exactly on the final ground state. In contrast, we show that if dephasing 
is added to the dynamics then there is a unique parametrization which maximizes the fidelity of 
the final state with the target ground state. The optimizing parametrization solves a variational 
problem of Lagrangian type and has constant tunneling rate along the path irrespective of the gap. 
Application to quantum search algorithms recovers the Grover result for appropriate scaling of the 
dephasing with the size of the data base. Lindbladians that describe wide open systems require 
special care since they may mask hidden resources that enable beating the Grover bound. 
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In the theory of adiabatic quantum control^ and quan- 
tum computation^, one is interested in reaching a target 
state from a (different) initial state with high fidelity, as 
quickly as possible, subject to given cap on the available 
energy. The initial state is assumed to be the ground 
state of a given Hamiltonian H{) and the target state 
is the ground state of a known Hamiltonian Hi. The 
two are connected by a smooth interpolating path in the 
space of Hamiltonians. The interpolation is denoted by 
Hq with q G [0, 1]. An example is the linear interpolation 

Hq = [1 - q)Ha + qHi, (0<g<l). (1) 

However, any smooth interpolation^ will do. We are in- 
terested in the optimal parametrization of the interpo- 
lating path. That is, a time-table for the path which op- 
timizes the fidelity of the time-evolving state at the end 
time with the ground state of the target Hamiltonian. 
We will show that when the evolution is unitary, there 
is no unique optimizer-there are plenty of them. How- 
ever, when dephasing is added to the dynamics there is 
a unique optimizer which we characterize and discuss. 

For the sake of simplicity we assume that the Hilbcrt 
space has a dimension N (finite) and that Hg is a self- 
adjoint matrix-valued function of q with ordered simple 
eigenvalues ea{q), so that 

N-l 

Hg=J2ea{q)Pa{q). (2) 
a=0 

Pa{q) are the corresponding spectral projections. 

A slow change of q tends to maintain the system in its 
ground state up to an error due to tunneling. We are in- 
terested in getting as close as possible to the target state 
within the time T allotted to traversing the path. The 
controls at our disposal arc a. The total time T and b. 
The parametrization of the path q{s) = (/^(s), s G [0, 1] 
for given e. Here s ~ et is the slow time parametrization 
and e = 1/T the adiabaticity parameter. 

The cost function is the tunneling Tq^^{l) at the end 



point, where Tq.£(s) is defined by 

T,,,(5) = l-tr(Po(gKs(s)). (3) 

Pq^e{s) is the quantum state at slow-time s which has 
evolved from the initial condition ^^^^(O) = Po(0)- 

A related but different optimization problem com- 
monly considered in quantum information is to optimize 
upper hounds on the tunneling^. The difference is that 
the cost function is evaluated not for a fixed, given in- 
terpolation, but for the worst case for any (smooth) in- 
terpolation between any two Hamiltonians belonging to 
certain classes. 

We consider two types of evolutions: (a) Unitary evo- 
lutions generated by Hq. (b) Non-unitary evolutions gen- 
erated by appropriate Lindblad generators L^. Since (a) 
is a special case of (b), both are of the form 

ep = Lq{p), (4) 

where ' = d/ds and 

M 

L{p) = -i[H, p] + ^ (2r,pr* - r;r,p - pr;r,) (5) 

3 = 1 

with Tj, a-priori, arbitrary. Adiabatic evolutions are 
a singular limit of the evolution equations since e hits 
the leading derivative. Unitary evolutions are generated 
when Fj = 0. 

In the case of unitary evolution the optimization prob- 
lem has no unique solution, on the contrary, optimizers 
arc ubiquitous. More precisely: 

Theorem 1 Let 

2Hq = g{q) ■ a (6) 

he any smooth interpolation of a 2-level system where a 
is the vector of Pauli matrices and g{q) a smooth, vec- 
tor valued function with a gap, \g{q)\ > go > 0; let 
£ / go he small. Then, in a neighborhood of order e of 



2 



any smooth parametrization, there are many non-smooth 
parametrizations with zero tunneling and therefore many 
smooth parametrizations with arbitrarily small tunneling. 

We shall sketch the main idea behind the proof. 
Consider a discretization of any given parametrization 
to (slow) time intervals of size 27:6 /go. In each in- 
terval one can find a point q* , such that the time- 
independent Hamiltonian Hq* acting for appropriate 
time T < 27r/|g(g*)| < 27r/5o, will map the image on 
the Bloch sphere of the starting point q- to the image of 
the end point q^. This says that there are many (non- 
smooth) paths, labelled by the continuous parameter sq 
in Fig. [U that map the instantaneous state at the ini- 
tial end point to the corresponding state at the final end 
point. These paths have zero tunneling. The existence 
of q* follows from the geometric construction in Fig. [1] 
g((7*) is a point of intersection of the path with the equa- 
torial plane orthogonal to g((7+) — g(<;-)- The resulting 
parametrization differs from the original one by at most 
(sup^ |9(s)|) • 27r£/(7oi as seen from the mean- value theo- 
rem. 



= Yl,a=Q \/7ja-^'" "^'hcrc ^ is a rectangular, MxN, 
matrix (without loss, M ~ iV^ — !)• It follows that de- 
phasing Lindbladians have the forro^: 

C{p)=~l[H,p]+Y,2^baPapPb-Y.^aa{Pa.p]. (7) 
a, 6 a 

where < 7 is a positive matrix. Time-dependent de- 
phasing Lindblad operators^ are then defined by setting 
H Hq and Pa Pa{q) and 7 -f{q) . 

The motion of ker Cq with q can be interpreted geomet- 
rically as follows: The space of (unnormalized) states is 
a fixed N'^ dimensional convex cone. The normalized in- 
stantaneous stationary states are a simplex whose ver- 
tices are the instantaneous spectral projections Paiq)- 
This simplex rotates with q like a rigid body, since the 
vertices remain orthonormal, tr(PaPb) = Sat and the mo- 
tion is purely orthogonal to the kernel, tr(P^Pb) = 
where P^ = dPa/dq. This follows from the fact that for 
orthogonal projections P^ is off-diagonal 

PLiq)-Y.P^(l^Pai'l)Pciq)- (8) 

b^c 




ex 



FIG. 1: Left: g± are the images on the Bloch sphere of the 
end points of an interval of size 0{e) of a given parametriza- 
tion (blue). The intersection of the associated interpolating 
path with the equatorial plane (shaded) determines the point 
q* and thereby the axis of precession g(q*) (red) that maps 
the instantaneous state at the initial end point to the corre- 
sponding state at the final end point. Right: A non-smooth 
interpolating path that takes the instantaneous eigenstate at 
the beginning of the interval to the instantaneous eigenstate 
at the end of the interval with no tunneling. 

Dephasing Lindblad operators belong to a special class 
of Lindblad operators which share with unitary evolu- 
tions the existence of N stationary states. (In contrast 
with generic Lindblad operators that have a unique equi- 
librium state.) More precisely, £ is a dephasing Lind- 
blad operator, if all the spectral projections Pa of H 
are stationary states, namely Pa € ker£. This is the 
case when [rj,iJ] = 0, and the condition is also nec- 
essary when H has simple eigenvalues, as can be seen 
by expanding ti{PaC{Pa)) — 0. We can then write 



An adiabatic theorem for dephasing Lindblad opera- 
tors can be inferred from^. It says that the solution p^^ 
of the adiabatic evolution, Eq. Q, for the parametriza- 
tion q{s) and initial condition p^°i(0) = Pa{0), adheres 
to the instantaneous spectral projectior^i^ 



p'q':Us)=Pa{s)+Oie), (s>0). 



(9) 



For the sake of writing simple formulas we shall, from 
now on, restrict ourselves to the special case where the 
positive matrix j{q) > of Eq. ([7]) is a multiple of the 
identity 

Cq{p) = ~t[Hq,p]^j{q)Y,P,{q)pPk{q)- (10) 

Our main results follow from: 

Theorem 2 Let Cq be the dephasing Lindblad of 
Eq. nop , and pq^^ a solution of with initial condi- 
tion p{Q) = Po(0) for the parametrization q{s). Assume 
a gap condition ea{q) 7^ eb{q), {a ^ b). Then the tunnel- 
ing defined by Eq. (0), is given by 



r,,,(l) = 2e / M{q)q^ds + 0{e^), 
Jo 



where the q dependent mass term 



7(g)tr(P,P^^) 



(eo(q)-e,(g))2+ 72(g) 



> 



(11) 



(12) 



is independent of the parametrization. Po(?) denotes a 
derivative with respect to q and q{s) with respect to s. 
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In the special case of a 2-level system, Eq. (|6]), where g{q) 
is a 3-vector valued function parametrized by its length 
dg{q) ■ dg{q) ~ {dqY the "mass" term of Eq. (|12|) takes 
the simple form 



Af(g) 



7(9) |gT(g) 



4 9^{q)+l'-{q) 



(13) 



|g'| is the velocity w.r.t. q on the Bloch sphere ball and 
5(9) = |g(9)l is the gap. 

Remark: For a 2-level system undergoing unitary evo- 
lution a similar variational principle to Eq. (jlip . but with 
a different M[q), was proposed, as an ansatz, in^i for the 
purpose of determining an optimal path, rather than an 
optimal parametrization of a given path. 

Before proving the theorem let us discuss some of 
its consequences: Note first, that the tunneling rate, 
2eM{q)q^ > 0, is local and uni-directional. It follows 
that whatever has tunneled can not be recovered, in con- 
trast with unitary evolutions. Eq. (fTT|) has the standard 
form of variational Euler-Lagrange problems with a La- 
grangian that is proportional to the adiabaticity e and 
with the interpretation of kinetic energy with position 
dependent mass. This variational problem has a unique 
minimizer qoi^s) in the adiabatic limit, in contrast with 
the case for unitary evolutions, which by Theorem [1] has 
no unique minimizer. 

Since the Lagrangian is s independent qo{s) conserves 
"energy" and the tunneling rate is constant along the 
minimizing orbit. This gives a local algorithm for opti- 
mizing the parametrization: Adjust the speed g(s) to 
keep the tunneling rate constant. The optimal speed 
along the path is then 



q 



M{q) 



(14) 



where r > is a normalization constant. This formula 
quantifies the intuition that the optimal velocity is large 
when the gap is large and the projection on the instan- 
taneous ground state changes slowly. The optimal tun- 
neling, Tmin, is then 



r„,in-2er + 0(e2), = / dq^M{q). 

Jo 



(15) 



This formulas will play a role in our analysis of Grover 
search algorithm. 

We now turn to proving Theorem [2] Evidently 



l-tr(Pop,,,)(l) 



d 



Now, the identity, 

C*iPaAPb) = {iiea-eb)-j)PaAPb, (a^b) (18) 
together with Eq. shows that 

PaPiPb 



«(ea - Cbj - 7 



(19) 



solves the equation 



P^{q)=C;{Xiq)) (20) 
Substituting this in Eq. p7|) gives the identity 



ds 



ty{Po{q)pqAs)) - etr(X(g) p,,,(s)) 9(s). (21) 



Integrating by parts the last identity gives an expression 
involving p but no p. This allows us to use the adiabatic 
theorem and replace p by P + 0(e). We then undo the 
integration by parts to get Theorem [2l 

In the theory of Lindblad operators H and Tj of Eq. ([5|) 
can be chosen independently. However, as we shall now 
show, if one makes some natural assumptions about the 
bath, the dephasing rate 7 of Eq. ([rU]) is constrained by 
the gaps of H. 

To see this we turn to quantum search with 
dephasingiiiS. Grover has shownii that 0{\/N) queries 
of an oracle suffice to search an unstructured data base 
of size ^ 1. The adiabatic formulation of the problem 
leads to the study of a 2-level system with a small gap 
given by^ii^ 



Aq)-^^-^ + il-2qf 

and large velocity on the Bloch sphere 



(22) 



(23) 



The time scale r, which determines the optimal tun- 
neling, can be estimated by evaluating the integrand in 
Eq. psp at its maximum, q = 1/2, and taking the width 
to be 1/^/N. This gives 



= O 



(24) 



to leading order in the adiabatic approximation. 

The adiabatic formulation^ fixes the scaling of the min- 



— tr(Po((7)pq^e(s)) ds . (16) imal gap go ^ but does not fix the scaling of the 



Using Eq. the defining property of dephasing Lind- 
bladians, Cq{PQ{q)) — 0, and by Eq. ([7]), the concomitant 
'C*(Po('7)) = 0, one finds 

^tr{Po{q)p,As)) ^ tr(P^(g)p,,,(s)) q{s) . (17) 



N 

dephasing rate 7 with N . We shall now address the issue 
of what physical principles determines the scaling of the 
dephasing with N . To this end we consider various cases. 

The regime 7 ^ e is outside the framework of the adi- 
abatic theory described here, but is close to the unitary 
scenario;^. For the adiabatic expansion and Eq. to 
hold e ^ 7. This means that in case of small dephasing. 
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7 ^ go, the allotted time, T » 7~"^ » 0{^/N), is longer 
than Grover search time. For such times the theory de- 
veloped here can be used to estimate the tunneling, but 
it is not appropriate for optimizing the search time. To 
optimize the search time one needs to study bounds on 
the tunneling rather than a first order term in e. 

When dephasing is comparable to the gap, 7 ~ 170, one 
finds M(l/2) ^ and from Eqs. ((Ml [22) one recovers 
Grover's result for the search time 

r^O^-^)=OiVN). (25) 

Finally, consider the dominant dephasing case: 7 3> .go- 
Here M ^ / gl and from Eqs. ([241 [22]) one finds 

r = 0(7-'). (26) 

If 7 scaled like 7 - iV-"/^ 1 > a., then T =0 (iV"/^) 
which seems to beat Grover time. 

The accelerated search enabled by strong dephasing 
is in apparent conflict with the optimality of Grover 
boundi^ii^: Consider the Hamiltonian dynamics of the 
joint system and bath, which underlies the Lindblad evo- 
lution. By an argument ofi^ for a universal bath, the 
Grover search time is optimal. How can one reconcile 
Eq. ([26]) with this result? Before doing so, however, we 
want to point out that Eq. ([^5]) is not an artefact of 
perturbation theory: While Tmin = 2£t is valid in first 
order in e, an estimate Tmin < er, with r as in Eq. (|T5|). 
remains true for all e provided 7 ^ 5o- 

The resolution is that a Markovian bath with 7 3> go 
can not be universal and must be system specific: The 



bath has a premonition of what the solution to the prob- 
lem is. (Formally, this "knowledge" is refiectcd in the 
dephasing in the instantaneous eigenstates of Hq.) Lind- 
bladians with dephasing rates that dominate the gaps 
mask resources hidden in the bath. This can also be seen 
by the following argument: Dephasing can be interpreted 
as the monitoring of the observable Hq. The time-energy 
uncertainty principloi^ says that if Hq is unknown, then 
the rate of monitoring is bounded by the gap. The accel- 
erated search occurs when monitoring rate exceeds this 
bound, which is only possible if the bath already "knows" 
what Hq is. When Hq is known, the bath can freeze 
the system in the instantaneous ground state arbitrarily 
fast. Consequently, the Zeno effecli^ then allows for the 
speedup of the evolution without paying a large price in 
tunneling. 

The formal theory of Lindblad operators allows one to 
choose the operators, H and Tj in Eq. ((51), independently. 
From the discussion in the last paragraph one learns that 
one must exercise care in using Lindbladians for systems 
that are wide open. Markovian baths which are univer- 
sal, i.e. oblivious of the state of the system, give rise to 
dephasing Lindbladians, with dephasing rates that are 
bounded by the spectral gaps of the system. 
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